arXiv: 1505.07226v2 [astro-ph.CO] 23 Sep 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 24 September 2015 (MN style file v2.2) 


On the Physical Requirements for a Pre-Reionization Origin of the 
Unresolved Near-Infrared Background 

K. Helgason*^, M. Ricotti^, A. Kashlinsky^’'^, V. Bromm^ 

^Max Planck Institute for Astrophysics, Karl-Schwarzschild-Str. 1, 85748 Garching, Germany 
^Department of Astronomy, University of Maryland, College Park, MD 20742 

^ Observational Cosmology Laboratory, Code 665, NASA Goddard Space Flight Center, Greenbelt MD 20771 
^SSAI, Lanham, MD 20706 

^ Department of Astronomy, University of Texas, Austin, TX 78712 


24 September 2015 


ABSTRACT 

The study of the Cosmic Near-Infrared Background (CIB) light after subtraction of re¬ 
solved sources can push the limits of current observations and yield information on galaxies 
and quasars in the early universe. Spatial fluctuations of the CIB exhibit a clustering excess 
at angular scales ~ 1 ° whose origin has not been conclusively identified, but disentangling 
the relative contribution from low- and high-redshift sources is not trivial. We explore the 
likelihood that this signal is dominated by emission from galaxies and accreting black holes 
in the early Universe. We find that, the measured fluctuation signal is too large to be pro¬ 
duced by galaxies at redshifts z>8, which only contribute ^0.01—0.05 nWm“^sr“^ to the 
CIB. Additionally, if the first small mass galaxies have a normal IMF, the light of their age¬ 
ing stars (fossils) integrated over cosmic time contributes a comparable amount to the CIB as 
their pre-reionization progenitors. In order to produce the observed level of CIB fluctuation 
without violating constraints from galaxy counts and the electron optical depth of the IGM, 
minihalos at z>12 must form preferably top-heavy stars with efficiency /*> 0.1 and at the 
same time maintain a very low escape fraction of ionizing radiation, /esc<0.1%. If instead 
the CIB fluctuations are produced by high-z black holes, one requires vigorous accretion in 
the early universe reaching pacc^lO®M 0 Mpc“^ by z~10. This growth must stop by z^6 and 
be significantly obscured not to overproduce the soft cosmic X-ray background (CXB) and its 
observed coherence with the CIB. We therefore And the range of suitable high-z explanations 
to be narrow, but could possibly be widened by including additional physics and evolution at 
those epochs. 
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1 INTRODUCTION 

The Cosmic Near Infrared Background (CIB) contains radia¬ 
tion that has been built up throughout the cosmic history, in¬ 
cluding highly redshift ed emission from the pre-reionization era 
(see e.g. review by iKashlins^ l2005h . Measured brightness 
fluctuations in deep source-subtracted CIB maps have estab- 
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iMatsumoto et alj 1201 ll : iKashlinskv et all l2012h . Its extragalactic 
nature is further supported by its isotropy (now mea sured in several 
different parts of the sky, IKashlinskv et alJ ll2012h l: 2) the signal 
is also inconsistent with the contribution from kno wn galaxy pop¬ 
ulatio ns at z<6 extrapolated to faint luminosities jHelgason et alJ 
I2OI2I) ; 3) the amplitude of the fluctuations i ncreases towards 
shorter wavelengths showing a blue spectrum jMatsumoto et alJ 
I2OIII : ISeo et alj|2015h : 4) 4) there is no evidence yet for the fluc¬ 
tuations correlating significantly with the mask or the outer p arts 
of removed sources jKashlinskv et al.l2005l : lArendt et alj2010h : 5) 
the large-scale clustering component does not yet appear to start 
decreasing as the small-scale shot noise power is lowered; 6) the 
clustering component of the fluctuations shows n o correlation with 
faint H ST/ACS source maps down to mAB~28 jKashlinskv et alJ 
l2007ch . 7) the unresolved CIB fluctuations at 3.6/im are coherent 
with the Cosmic Far-Infrared Background (200,350,500^m) which 
can largely be expl ained by unresolved galaxies at low-a and their 
extended emission dXhacker et al.ll20f^ . 8) the source-subtracted 
CIB fluctuations are coherent with the unresolved soft Cosmic 
X-ray Background (not detected in harder X-ray bands of >2keV) 
j^Dpelluti et all [201311 which can partly be ac counted for by 
unres olved AGN and X-ray gal axies at low-z 1 Helgason et alJ 
l2014h . A recent measurement by IZemcov et al. |j20li finds the 
blue spectrum continuing to shorter wavelengths but appears to 
be in conflict with earlier //SF/NICMOS and 2MASS fluctuation 
studie s at the same wave lengths from 2MASS Kashlinsky et alJ 
l2002h and //ST/NICMOS jThompson et al.ll2007alk 

It was suggested that the era of the first stars could have 
left a measurable imprint in the CIB, b oth its mean level 
(ISantos_et_alj_ 200^ and its anisotropies dCoorav et alJ 1 20041 ; 
IKashlinskv et alJ 2004ll . This was followed by more detailed 
studie s on the nature of these populations jSalvaterra et"^ 
20061; IS alvaterra & Ferrar j l2006l; iF ernande z & Komatsul l20od ; 


Kashlinsky et al. 2007bl ; Fernandez et al. 2010h . This was moti¬ 
vated by the expectation that the first objects were i) individually 
bright with a short epoch of energy release, ii) highly biased as 
they form out of rare density peaks, and iii) radiate strongly in 
UV/blue being red shifted into to day’s near-IR part of the spec¬ 
trum (reviewed in iBrom i^ [20l3h . Despite being consistent with 
all observed properties however, the high-2 origin of the source- 
subtracted CIB signal continues to be debated. This is in part be¬ 
cause of the lack of a robust redshift determination of the signal and 
because recent models of early galaxy popul ations 2>6 have failed 
to produce sufficien t CIB fluctuation p ower jFemandez et alj2012l; 
ICoorav et alJl20123 ; lYue et ^l2013al) . This has motivated alterna¬ 
tive hypotheses for t heir origin, such as in a diffuse intrahalo lig ht at 
low/intermediate-2 jCoorav et alj|2012bl ; IZemcov et alj|2014l) . At 
the same time, it was suggested that accretion by direct collapse 
black holes (DCBFI) at high-2 can provide an explanation w hich 
also fully accounts for the CIBxCXB signal jYue et alj2013bh . 

We evaluate the physical requirements for a pre-reionization 
origin of the CIB fluctuations based on the latest observational 
insights. We carefully quantify the clustering excess in Section 
In Section we use an analytic approach to model the CIB 
contributions from early stellar populations, both metal-free and 
metal-enriched. In Section |4] we consider accreting black holes 
and associated gaseous emission. We derive the required star for¬ 
mation and accretion rates and discuss their implications for the 
CXB. In this paper we use the standard ACDM cosmology with 
parameters {h, Dm, Da, Ws , cr8)=(0.678,0.308,0.692,0.968,0.829) 
jPlanck C ollaboration et al. ll2013l) . All magnitudes are in the AB 
system i Oke & Gunnlll983n. 


2 QUANTIFYING THE CIB FLUCTUATION EXCESS 

Brightness fluctuations in the CIB can be written 5T’(x)=T'(x) — 
(F) where T’(x) is the sky brightness at the 2-D coordinate x and 
(F) is the mean isotropic flux. We describe the fluctuation field in 
terms of the power spectrum as a function of the angular wavenum¬ 
ber q, defined as P(q)=(|(5q|^), where 5q is the 2-D Fourier trans¬ 
form of the CIB fluctuation field, 5F{'x.). On the angular scale 
9=2111q, the root-mean-square fluctuation in the CIB can be writ¬ 
ten as, (5Te = (5F^(x))^^^~[g^P(g)/27r]^/^. The fluctuations are 
therefore determined by the flux of the underlying sources and how 
they cluster on the sky 

5Fe=FciB^e (1) 

where the information on the clustering as a function of angular 
scale is contained within the fractional fluctuations Ag. At high- 
2, the sources trace increasingly rare peaks in the density field and 
Aci can become ~0.1 at arcminute scales. Current fluctuation mea¬ 
surements exhibit a signal above the noise which flattens to an ap¬ 
proximately constant rms value towards large angular scales (see 
Figure [TJ. In this paper we refer to the clustering as the average of 
this large scale value Aci=(A(3'—30')). When defining our CIB 
fluctuation excess, we focus on NIR fluctuations at 2<A<5/rm 
from AKARI/IKC and SpizlerflRAC. This is where current mea¬ 
surements show a mutually consistent large scale signal which at 
the same time are deep enough such that the clustering component 
is not sensitive to further decreasing the shot noise from remaining 
galaxies. 

In order to get a single representative value for the mea¬ 
sured large scale CIB fluctuations from high-2, we consider a sim¬ 
ple three-parameter model consisting of the sum of fluctuations 
from unresolv ed galaxy population s at the measurement detec- 
tion threshold jHelgason et alj2012h and ACDM power spectrurrQ 
jEisenstein & HrJI 19981) . projected to 2=10: 

\ \ 

+c.5F,1,(A). ( 2) 

We explore the parameter space by least-squares fitting the avail¬ 
able data from Akari and Spitzer (at l'—50') in 10® Markov-Chain 
Monte Carlo steps. We use Metropolis-Hasting acceptance with 
chain bum-in of 10® steps and priors on the parameters — l<loga< 
1, —3<logb<3, —2<logc<5. As each measurement has different 
source detection threshold we allow the last term in the above Equa¬ 
tion to vary at each measurement depth such that a is a free pa¬ 
rameter in six different measurements (making the effective num¬ 
ber of parameters eight). For the first term, we obtain a=1.79lg’4f 
and 6=1.77ltQ'5®. We display the 68% confidence regions of the 
best-fit in Figure [T] for which we obtain x^/dof=1.62. The qual- 
ity of the fit is limit ed by t he slight disagreement i n the data of 
ICoorav et alj j2012bh versus IKashlinskv et alj j2012h . particularly 
at 4.5/im where our single power-law spectral slope favors the lat¬ 
ter. It is worth noting that the best-fit excess fluctuation signal, 
does not exhibit a Rayleigh-Jeans type spectrum, ~A“® 
as indicated by measurements of the net large-scale signal. The 
reason is the contribution of remaining low-2 galaxies: 1) because 
the abundance of low-2 galaxies is greater at the typical detection 
thresholds at <2fim i.e. their number counts are steeper, and 2) be¬ 
cause experiments at smaller wavelengths happen to be shallower 
than those at >2fim such that less of the low-2 galaxy contribution 

^ The power spectrum is normalized at the 10' peak such that i5T3.6(jm = 
0.043nWm~^sr“^. 
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Figur e 1. The latest measurements of the mean-square CIB fluctuations in the 2-5trm range. AKA RI/IRC measurement s are shown as red circles JSeo et alJ 
I 2 OI . 5 I 1 whereas Spitzer/IRAC are shown in blue circles iKashlinskv et al.ll2012ll and green squares iCoorav et alJl20123l . The solid c urve shows the best-fi t 
ACDM power spectmm projected from 2=50 to 2=10 whereas the dashed lines correspond to the unresolved galaxy populations from iHelgason et S] i2012ll 
fitted to each individual measurement (see text). The fit gives x^/dof=l-62. The upper horizontal axis shows the angular scale in arcseconds. The large scale 

averaged fluctuations in each band are 5F((3'—30'))=(0.11,0.068,0.056,0.037)^|p'°g^’p °j 2 ’o o 06 'o 009 ) tiWm“^sr“lat (2.4,3.2,3.6,4.5)/im. 


is removed. The SED of the best-fit excess can be seen in Figure|^ 
Integrating the first term in Equation over 2-5/im we obtain the 
excess fluctuation 


3 STELLAR SOURCES 
3.1 High -2 Galaxy Populations 


5Fcm = 



&F'^ 


dA 

A 


1/2 

—0 070 + 0-023 

—U.U/Z_o.o20 


nW m ^sr ( 3 ) 


This c ompares well with the analytical estimate in iKashlinskv et all 
) l2015h . We will refer to this value of (5Tcib= 0.072 nWm“^sr“^ 
as our requirement for any high -2 model to reproduce the data in 
Figure [T] We also make our equations directly scalable with Ad 
such that the associated CIB flux TbiB=5TciBcan be easily 
compared with measurements of the isotropic CIB, both direct and 
those derived from TeV blazars. 

At a given near-IR wavelength, the isotropic CIB flux is re¬ 
lated to the comoving specific emissivity jv ( 2 ) per unit volume of 
the sources 


where v—VaheU-Vz) is the rest frame frequency. If the CIB was 
released during the first Af~ 500 Myrs, it follows from this Equa¬ 
tion that the luminosity density pi^=vji, at optical wavelengths 
0-2<Aci b/( l+,2)5-,0-5ftm must reach a representative value of 
(see also iKashlinskv et alj2007bh 


Pl=1.5x 10®I/Q Mpc ^ 




Af 


500 Myr 


( 5 ) 

where we have substituted Ad=(5TciB/TciB- This luminosity 
density is notably higher than the output during the peak of star 
formation history at 2~1—2. If the CIB fluctuations originate at 
high- 2 , the luminosity density must have been substantially higher 
in the early universe regardless of the nature of the sources. In the 
following Sections we investigate the basic astrophysical require¬ 
ments for both stellar and accretion powered emission. 


An increasing number of galaxies are being detected out to 2~10 
as the deep Hubble pro gram are being pushed to the limits of the 
instrument capab i lities i Bpuwenset^ 2007 : McLure et al.l 200^ : 
lEllis et all I 2 OI 3 I : iFinkelstein et al. 2014 : Bouwens et al.UOlA ). 


While these galaxies compose merely the tip-of-the-iceberg at these 
redshifts, it is of interest to estimate the expected contributions of 
the entire population to the unresolved CIB. Here, we present a 
simple forward evolution model of the conditional luminosity func¬ 
tion which we tune to fit observations at 2^8—10. We assume that 
galaxies form in halos with a SFR proportional to their collapse rate 
at 2 


p*(M,2)=/*^+Mn(>M^in,f) (6) 

i LM CLt, 

where /* is the average fraction of baryons in collapsed halos that 
are processed into stars an d n(M,t) is the evo lving halo mass func¬ 
tion, which we adopt from iTinker et al.l ( l2008h . In order for star for¬ 
mation to take place, halos must reach a sufficient size, Mmin, to 
shock heat the baryons to the virial temperate allowing for efficient 
gas cooling. The relation between the mass of a halo and its virial 

temperature is Minin=3x 10®^ where 

we have assumed a mean molecular weight of 1.22 for neutral pri¬ 
mordial gas composed of H and He. For the discussion of high- 
2 galaxies we will adopt the limit of Tvir=40,000iT, somewhat 
higher than the threshold for cooling via atomic hydrogen. Star for¬ 
mation in smaller halos will be addressed in the following subsec¬ 
tion. Luminosity and spectra are assigned to star forming galaxies 
using the Yggdrasil model, a population synthesis code designed to 
mode l high -2 systems contai ning varying mixtures of PopII and III 
stars llZackrisson et al J 2 OI ih . This model includes the nebular con¬ 
tribution from photoionized gas and extincti on due to dust, w i th sin- 
gle age stel lar popula t ion (S SP) taken from iLeitherer et ^ ) ll999h 
(PopII) and lSchaererl (l2002l) (PopIII, with no mass loss). Because 
observed Lyman break galaxies show no evidence of anything other 
than metal-enriched star formation, we adopt a universal Kroupa 
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Figure 2. The single age stellai' population spectra from the Yggdrasil code. 
The three rows show energy spectra for different metal-free/poor IMFs. 
From top to bottom panels: metal-free '^500 Mq stars radiating close to 
the Eddington luminosity, heavy IMF with characteristic mass ~100Mq, 
intermediate IMF '^IOMq and metal-poor stars with a normal Kroupa IMF 
(see text for details). The left and right panels show the SSP spectra for the 
two limiting cases of /esc= 0 and 1 respectively. The three curves in each 
of the panels shows the single age population at 0 and 3.6 Myr (top two) 
and 0, 10, 100 Myr (bottom two). 


IMF (0.1—IOOMq) with a metallicity of Z=0.0004, characteris¬ 
tic of PopII stars (we call this IMFi in the discussion below and 
Figure[2j- The volume emissivity of halos with masses between M 
and M-\-dM at 2 : is obtained by convolving the SSP with the star 
formation occurring prior to 2 : 

/ OO I / 

(7) 

where jCv{tage) is the aging spectral template shown in Figure |2] 
At any given epoch, each halo includes the instantaneous emission 
from newly formed stars as well as older populations from earlier 
episodes of star formation. The luminosity function is d?(L)dL= 
<E>(M)dM where the relation between mass and luminosity is 
L{M)=i/ju{M,z) / {Mn{M,z)). The CIB flux production history 
seen at frequency izohs as a function of halo mass is 

(8) 


where is evaluated at j/=i.'obs(l+2)- The net CIB is then simply 
Tcib= // f{M,z)dMdz. This can be used to derive the angular 
power spe ctrum of CIB fluctuations via projection of the source 
clustering l lLimbeilll953h 


Pi<})= II f{M,z)b{M,z)dM 


Paiqd ^,z)dz (9) 


Figure 3. Upper: The SFR in halos (Tvir>4xlO^K, gr een) and miniha¬ 
los (T vir>10^K, blue). The SFR data points, compiled in iRobertson et alJ 
are derived from the UV luminosity function integrated down to 
O.OOIL*. The solid line shows the maximum likelihood SFR history from 
the same work. These connect well with our modeled SFR at higher- 2 : as¬ 
suming /*=0.005 Bottom: The build-up of large scale averaged CIB fluc¬ 
tuations from galaxies forming stars with IMFi and efficiency /*=0.005. 
The solid and dashed lines correspond to /esc= 0 and 1 respectively. Fluc¬ 
tuations from I0W-2: galaxy populations remaining after re moving galaxies 
brigh ter than r723.6^m=25 are shown as the black solid line IMelgason et alj 
l2Q12n . The reference level of the observed large scale fluctuations are shown 
as the hatched region (grey). 


where and d{z) is the comoving dis¬ 

tance. The clustering is described b oth in terms of the evolv ing mat¬ 
ter power spectrum in 3D, Ps{k,z) teisenstein & Hull998h and the 
mass dependent halo bias, b{M^z) fsheth et a i Booiir In this de¬ 
scription, the bias is coupled to the brightness distribution of halos, 
eliminating the need for assumptions on the halo occupation num¬ 
ber of galaxies. We also note, that variations in the cosmological 
parameters, in particular as, affect the amplitude of the power spec¬ 
trum and adds uncertainty to our modeling. Equation|9]is equivalent 
to the 2-halo t erm describing the corr elation between central halofl 
llSeliakl200d:ICoorav & Sh^l2002h . 

The large scale CIB fluctuations resulting from our modeled 
high -2 galaxies are shown in Figure We have chosen a constant 
star formation efficiency (in both M and z) in such a way to ob¬ 
tain a good agreement with the observed star formation rate and 
UV LF at high- 2 . This value we find to be /*=0.005. Figure |4] 
compares the relative abundance of high -2 galaxies (green) with 
the source density of all z<8 galaxies (blue). It is immediately 


^ In this paper we neglect the 1 -halo term from high -2 sources since it is 
always negligible compared to the 2-halo term at the relevant angular scales 

03'). 
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Figure 4. Deep source coun ts at 1.25Atm (left) and 1.6^™ (right). The blue data points represent galaxy counts, mostly from 3 at these deep magnitudes. 

The galaxy counts model of lHelgason et al.l i2012h is shown as blue solid line and extrapolated to faint magnitudes whereas the green line shows the counts 
for galaxies at high-z (/*=0.5%, IMFi), with the Lyman dropouts placing them at z>8 (left) and z>10 (right). The comparison of low-z and high-z 
counts suggests that the faint-end of the high-z galaxy population will only become significant at >32 mag. The dashed lines show the model that would 
required to reproduce the CIB fluctuations, (SFcib, at these redshifts. The fact that these models are in gross conflict with both the LF data and the deepest 
counts suggest that any such signal must either come from higher z or be confined in isolated minihalos at the faint-end. The inse ts show the luminosity 
function log(mag~^ Mpc~^) vs Myv of our default high-z galaxy model (/*=0.5%, IMFi, /esc=12%) tuned to fit measurements iFinkelstein et alJ2014l : 
iBouwens et alj20i3l . 


clear that compared to current measurements, star forming galax¬ 
ies with a normal IMF are underdominant at least out to 32 mag 
and, with these assumptions, do not reach the measured fluctua¬ 
tion levels (see Figure O. This is true regardless of whether the 
LF is cut off at Muv——14 mag or extrapolated to much fainter 
systems. Simply increasing /* further would overproduce the LF at 
2 ~ 8 —10 and eventuall y also the faint NIR number counts (this was 
already pointed out by ISalvaterra & Ferr^l 20061) . In fact, fluctu¬ 
ations from faint unresolved galaxies at low-z produce a larger sig- 
nal than high - z galax ies. We therefore conc lude, in agreement with 
ICoorav et alj i2012ah and lYue et al.l (1201 3al) . that the measured CIB 
fluctuation levels cannot be reproduced by high-z galaxies with rea¬ 
sonable extrapolations of their evolving luminosity function. 

3.2 Minihalos and PopIII 

There are two ways of having more light produced by stars without 
overproducing the observed LF at z~ 8 —10. First, that an epoch 
of more vigorous light production took place before z~ 8 — 10 ; in 
other words, the era of CIB production already ended before this 
time. Second, the sources are outside the sensitivity limit of the 
deepest surveys i.e. they are found in numerous isolated halos that 
are intrinsically fainter than Mab— —17 mag. This would be ex¬ 
hibited in a rise in the faint-end of the LF, either by greater star 
formation efficiency in smaller halos or by heavier IMF. How¬ 
ever, numerical studies tend to sho w the opposite, that th e effi¬ 
ciency decreases wi t h lower masses llRicotti et 'Zl l2002tffl[200i; 
iBehroozi et alj2013l ; IO’Shea et al.l2015l) . In this subsection, we re¬ 

lax the assumption of a constant /*=0.005 and Tvir^TOjODOK and 
explore whether greater values can be accommodated where ob¬ 
servational constraints are not yet available. In other words, what 
happens if we, in addition to high-z galaxies > 1 O®M 0 , include the 
contribution of the first stars forming out of pristine metal-free gas 
in minihalos > 10 ®Mq? 

The IMF of the first PopIII stars is highly uncertain but is ex¬ 


pected to be biased towards high masses (iBrom ml l2013l) . It is also 
unclear whether stars forming under 772 -governed cooling (some¬ 
times called PopIII. 1) differ substantially from those forming in a 
second episode governed by atomic H line emission (PopIII.2). We 
therefore consider four different IMFs of the Yggdrasil model for 
stars forming in halos down to the smallest minihalos equivalent to 
the H 2 cooling threshold 1,000K 

IMFi: standard Kroupa in the 0.1 —lOOM© range 
IMF 10 : log-normal with characteristic mass of lOM© and dis¬ 
persion of IM© in the 1—500M© range 

IMFioo: power-law in the 50—500M© range 

IMF500: all stars are 500M©, near-Eddington 

For details, we refer to IZackrisson et alj 1 I2OI ih and references 
therein. Spectral evolution templates for these IMFs are shown in 
Figure|2 

Figurej^shows the average star formation efficiency required 
to produce the measured (5FciB=0.072nWm“^sr~^ at any given 
epoch for the three IMFs considered. It is clear that under these 
assumptions the CIB fluctuations can only be reproduced if mini¬ 
halos ~ 10 ®M© are allowed to continuously form massive stars 
with a high efficiencies /*>0.1. Furthermore, deep NIR counts 
limit this possibility to the most heavy IMFs where much of the 
energ y output is reprocessed into nebular emission, /esc =0 (see 
also iFernandez et alj l2010h . This is assuming that stars form 
continuously without mechanical or radiative feedbac k which is 
known to have strong imp act on both /* and /esc (e.g. Ijeon et all 
|2014 IPawlik et~^ 120151). In Figure | 6 ] we display the full CIB 
model with a simple description of chemical feedback which mod¬ 
els the PopIII-Po pII transition in terms o f semi-analytic mod- 
elling of SN winds llFurlanetto & Loebl2005t) which we adopt from 
iGreif & BrommI ( l2006l) 

PPopIII—Ppris(z)/)* (10) 

where p* is the net SFR from Equation and Ppris is the frac- 
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Figure 5. Upper. The fraction of mass in collapsed stimctures as a function 
of redshift. Solid and dashed lines correspond to halos with T vir>l,000K 
(>10^ MiT)) and >40,000K (>10® Mq) respectively f or the [Tinker et al 
200^ ) mass function. Dotted lines show the classic IPress & Schechtei 
19741) prediction for comparison. The grey regions show the scatter as¬ 
sociated with the 1-cr uncertainties in the Planck cosmological parameters. 
Lower: The star formation efficiency /* required to produce the reference 
level of the CIB fluctuations, SFqib^ by a given redshift, ^lend- The curves 
assume the entire stellar' population forming with IMFi, IMFio, IMFioo. 
IMF 500 (red. green,blue,orange) in all halos. The solid lines show the case 
where minihalos Tvir>l,000K are included whereas the dashed lines in¬ 
clude Tvir>40,000K halos only. We have set fesc=0 in all cases and note 
energy requirements become even greater for /esc>0. The colored regions 
show the combination of /* and ^end that result in the overproduction of 
deep NIR counts data and/or LF data for the IMFs in the same color scheme 
(see Figure[4j. These regions are thus forbidden unless the bright-end of the 
LF is suppressed e.g. if /* much lower in high mass halos than in low mass 
ones. 


tion of collapsed objects that are still chemically pristine. We dis¬ 
play two cases of PopIII-PopII transition in Figure standard 
and delayed enrichment corresponding to Zi/ 2 =^^ and 15 respec¬ 
tively where we have defined Zi /2 as the redshift at which half of 
all star formation is in metal-enriched modJE The PopII assumes 
IMFi whereas we display the cases where PopIII takes on IMFio, 
IMFioo> and IMFsoo- Only the most extreme models (/*>0.1, 
/esc^l, IMFioo, 50 o) come close to the fluctuation excess shown 
in the right panel as a grey region. 

3.3 Ionizing Photons 

High star formation efficiencies coupled with a heavy IMF imply 
a vigorous production of ionizing photons (>13.6 eV). These pho- 

^ The enrichment histories are modified by varying the radiation loss pa¬ 
rameter (see iGreif & Bromml2QQd) . 


tons must either be absorbed locally by neutral gas within the halo, 
or escape and be absorbed by still neutral parts of the IGM. The 
injection rate of ionizing photons can be derived self-consistently 
from the SFR combined with the stellar population synthesis mod¬ 
els at z>10 

nion{z)= f f dMdu. (11) 

7>13.6eV/h7 

For the production of i onizing photons at late r times z<10, we 
adopt the prescription of iRobertson et all (1201 Sh where fiion is de¬ 
rived from fitting the measured star formation rate history (see Fig¬ 
ure O. We calculate the ionization fraction of the IGM, Xion, by 
solving 

^ _/escflion(2:) Xion /i-i\ 

Xion } T ~ (tZ) 

{J^h) free 

where nion(z) is the production rate of intrinsic ionizing photons 
per comoving volume and (nn) is the average comoving num¬ 
ber density of hydrogen. The recombination timescale is frec = 
[Casn_ff(l-|-y/4W)(l-|-2)^] ^ where C=(n^H)/( uh)^ is the 
clumping factor of ionized hy drogen, aB(20 ,000K) is the case- 
B recombination coefficient dhlumme il ll994h and X=0.76 and 
Y—l—X are hydrogen and helium abundances. This allows us to 
derive the Thompson optical depth to the electron scattering 

Te=caT{nH) j Xioii(2) ^1+^^ (l-|-2)^^dz (13) 

where we assume Helium is singly ionized at 2>3 (t;= 1) and 
fully ionized at later times (rj=2). This can be compared with 
the latest values inferred from CMB polarization measurement s 
from Planck re=0.066±0.012 dPlanck Collaboration et al.ll2015h . 
which were found to be so mewhat lower than e arlier results from 
WMAP Te=0.088±0.014 dHinshaw et al.l2013h . In Figure|7](left), 
we show the integrated optical depth to electron scattering, Xe, as a 
function of redshift for our preferred IMFi model of high-z galax¬ 
ies. Forcing the model to reproduce the latest Planck measurements 
requires /esc = 12% and 5% for minimum host halo mass corre¬ 
sponding to 40,000K (~1O®M0) and lO^K (^IO^Mq) respec¬ 
tively. For the somewhat higher WMAP-derived Xe these values 
are /esc=32% and 12% respectively. We caution that, since we did 
not include any feedback effects which tend to surpress the star 
formation efficiency at small host halo masses, the average escape 
fraction may be somewhat higher. 

For a given SFR and IMF, the injection rate of ionizing pho¬ 
tons into the IGM is simply proportional to /* and inversely propor¬ 
tional to /esc. Figure|7](right) shows the combination of /* and /esc 
that give a fixed Xe=0.066 for our four IMFs. The Figure shows 
that in order to maintain the required /* and at the same time avoid 
reionizing the universe too early (and overproducing Xe), one must 
maintain a very low average escape fraction /esc)S0-l% at 2 > 10 . It 
follows that the >99% of the ionizing photons need to be absorbed 
by the local gas which itself can only constitute (!—/*)< 90% of 
the available baryons in these small halos. Several studies modeling 
the propagation of ionization fronts within halos hosting the first 
galaxi es show that the escape fraction increases towards smaller 


mass llRicotti & Shulll|200C 

; Johnson et al. 

20091) ■ For star forma 

tion efficiencies /'*>10“®, 

Ferrara & Loeb 

(I 2 OI 3 I) find that ioniz 


ing photons escape very easily with /esc—1 across a wide range 
of halo mass. Only the most massive halos, >10® Mq are able to 
confine their UV photons effectively. The escape fraction is also 
expected to increase towards hig her redshifts making it hard to jus- 
tify large /* in high-z minihalos dKuhlen & Faucher-Giguerel2012l : 
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Figure 6. The spectrum of the CIB flux, FQiB=yIi^ (left) and large scale fluctuations, 5FQiB = yJq^P/2'K (right) for models of 2:>8 early stellar emission. 
Fluctuations are averaged in the angular range 3^—30^ The observationally motivated model of high -2 galaxies, with /*=0.005 and IMFi is shown for 
Jgg(,=0.12 (black solid line) and /esc=l (black dashed) for comparison. The colored set of lines correspond to models with /*=0.1, IMFi representing 
PopII (solid), and three cases of PopIII IMFioJMFiooJMF 500 (dot-dash, dashed, dotted). The amplitude of both flux and fluctuations is directly scalable 
with /*. Orange corresponds to a PopIII-PopII transition with Zi/ 2 =^^ according to Equation 1 101 whereas blue illustrate the case of a delayed enrichment 
with (see text). All star formation is ended at 2 = 8 . Models shown in blue and orange are for /esc=0 since /esc>0.1% already violates reionization 

constraints explained in Section[33] Left: The data points is CIB from integrated galaxy counts which is naturally much greater than the high -2 conti'ibution. 
Right: The grey shaded region shows the best-fit CIB fluctuation excess obtai ned in Section[2]after the contributio n of low -2 galaxie s has b een accounted for. 
The data in the r ight panel is large scale averaged fluctuat i on me asurements of IZemcov et alJ42014h (filled circles), iMatsumoto et al.l j^llh (open diamonds), 
ISeo et al' ] | 20 il (filled diamonds) and iKashlinskv et"^ j2012h (filled squares). We do not show 2MASS and //^TTNICMOS measurements as they do not 
reach sufficient angular scales for proper comparison at 3'—30^ 




Figure 7. Left: The Thomson electron scattering optical depth integrated over redshift for our standard mo del of galaxy formatio n at 2 > 10 . The 2<10 
conti'ibution shown in light red is derived from the measured star formation rate history shown in Figure[^(see iRobertson et al .120151) . The red lines show the 
escape fraction tuned to match rg from Planck (lower) and WMAP (upper) which are indicated as horizontal dotted lines with blue regions showing the Icr 
uncertainty. Right: The mean escape fraction versus star formation efficiency in all collapsed halos for a fixed value of re=0.066 assuming the four IMFs in 
the same color scheme as in Figure[5] Dashed corresponds to a cooling threshold of Tvir>l,000K ('^lO® Mq) and dot-dashed to Tvir>40,000K {'^10^Mq). 
In order to reproduce SFqib ^nd at the same time be consistent with reionization, the combination of (/esc,/*) must be concentrated in the lower right comer. 
The arrows denote the minimum required /* are taken from Figure[5](lower) at 2 = 10 . 
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iMitra et albolsh . There are recent indications from determinations 
of the high -2 LF that the star formation effic iency is actually in¬ 
creasing at early times llFinkelstein et al.l2015h . While it is not clear 
whether this trend extends to the lowest mass halos where cool¬ 
ing relies on H 2 , observations of the stellar mass function in local 
dwarf systems suggest that small atomic cooling ha los at high -2 
could indeed be very efficient sites of star formation dMadau et alj 
|2014) . 


3.4 Supernova contribution 

Th e contribution of supernov a to the CIB was briefly discussed 
by ICoorav & Yoshidj j2004) who argued that SNe would re¬ 
main subdominant with respect to the stellar contribution. The 
net energy radiated in a core-collapse SN is i5sN~10®^erg out 
of which only 0.1-1% emerges as electromagnetic radiation. For 
comparison, a star burning 10 % of its initial hydrogen emits i?* = 
1 O®^(M*/M 0 )erg over its lifetime, out of which only a fraction of 
this energy will end up in the CIB, depending on the spectrum and 
redshift. A crude calculation can be made based on our model by 
assuming that every >8Mq star explodes as a supernova for our 
three IMFs above. Even if every SNe inputs a generous 10“*^ergs 
into the CIB, the contribution is always <0.1% of the stellar main- 
sequence contribution. It is therefore safe to ignore SNe as a major 
CIB contributor even in the case of some unusu ally bright SN-types 
(PISN, hypemova) that have been proposed dBarkat et alJ l l967l : 
llwamoto et al.lll99^ . 



Figure 8. The integrated 2-5/rm flux production when star formation is 
turned off at 2=10, normalized to unity at 2=10 (arbitrary units). The red 
curve corresponds to a standard IMFi containing ageing low-mass stars 
in all Tvir>10®K (>10 ®Mq) halos. The green and blue lines show the 
heavier IMFio and IMFioo respectively where the CIB contribution drops 
quickly after star formation is shut off as heavier stars die. The curves have 
been normalized at 2=10 for easier comparison. 


3.5 Fossils of the first galaxies 

Around the time of reionization, it is expected only systems with 
sufficient mass (> 10 ® —IO^Mq) will be shielded from photoheat¬ 
ing and continue to form stars. The majority of the already col¬ 
lapsed low-mass systems can nevertheless continue to contribute 
to the CIB through ageing stellar populations that continue to 
radiate throughout cosmic time. Locally , these systems are re¬ 
ferred to as the fossils of the first galaxies dRicotti & Gnedinll2005l : 
iBovill & Ricottil[2009l) . These will not be as intrinsically bright as 
their high -2 progenitors but can in principle accumulate substan¬ 
tial CIB by emitting at much lower -2 and for much longer span 
of time. Our formalism outlined above accounts for the ageing of 
populations and allows us to estimate the fossil contribution at later 
times. Approximating reionization as a step function, we turn off 
the SFR at 2=10 and allow the galaxies to passively evolve to 2 = 0 . 
Not all fossil stars will contribute to the source-subtracted CIB as 
many systems will be incorporated into larger halos hosting star 
formi ng galaxies that are ma sked and subtracted in a fluctuation 
study teovill & Ricottil201 lal fQ). For simplicity we do not account 
for any subtraction of the fossil contribution and the calculation 
presented here includes all star-forming progenitors at 2 > 10 . This 
should therefore be taken as an upper limit to the source-subtracted 
CIB for pre-reionization fossils. Note, however that this simple pic¬ 
ture has recently been challenged by recent evidence showing that 
many dwarf systems contin ued to form stars even well after the 
universe was reionized (e.g. IWeisz et al.ll2014h . 

Figure [ 8 ] shows the emission history of star forming galaxies 
( 2 > 10 ) and their remnant fossils ( 2 < 10 ). Shortly after SFR is set 
to zero, the flux from IMFio and IMFioo drops quickly as more 
massive stars die out but for IMFi the contribution of the remaining 
low-mass stars gradually rises towards lower -2 due to the ( 1 - 1 - 2 )“^ 
factor in Equation [4] Interestingly, when integrated over redshift 


the CIB contribution from fossils is ~70% of the net contribution 
whereas high -2 galaxies make up only ~30%. This suggests that 
the fossil contribution should be included in the net energetics of 
background radiation. However, considering that the low -2 part is 
an upper limit for the fossil contribution, this is still insufficient to 
account for the measured JFcib. It will only become important if 
the corresponding high -2 contribution is shown to be significant. 


4 ACCRETION POWERED EMISSION 

CIB originating from accretion onto black holes offers a more en¬ 
ergetically favorable scenario due to the more efficient conversion 
of mass to energy. Such a scenario is particularly attractive for 
explai ning the CIB in light of the observed coherence with the 
CXB dCappelluti et al.ll2013h and the expectation that the progen¬ 
itors of supermassive black holes (SMBH ) in the local universe 
were growing rapidly during this epoch (see IVolonteri & BellovarvI 
I2OI2L for review). With an average radiative efficiency of ac¬ 
cretion t=L/M(?, a black hole will increase its mass by M= 
(1—e(M,f))Min where Min is the net incoming mass accretion 
rate. Throughout this section we will assume a constant efficiency, 
e(M,f)=const with a caveat that relaxing this assumption may al¬ 
ter the numerology that follows. The net energy density radiated 
by a black hole pop ulation can be related to mass density accreted 
before redshift 2end ( ISoltanlll982h 

l_e /■°° dt f 

Pacc{Zend) = ^ / Pju{z)du (14) 

which is independent of the individual masses of the black hole 
population. The energy output is related to the accretion reate via 
pL=ef)a,cc(? for which the CIB can be written (see Equation|4ll 

-Fcib = — / epaccrfi / -r— (15) 
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where bi, is the average emerging spectrum, normalized such that 
fb^diy—1 and taken to he independent of mass and time. When the 
last term is roughly constant or slowly varying at early times (in 
particular for a blue spectrum sloping up with shorter wavelengths) 
the net CIB from accretion processes is simply proportional to the 
net mass accreted by redshift jZend i-e- ^ciBOcpacc- In other words, 
with these assumptions the CIB is mostly independent of both the 
growth history and the individual masses of the BHs. It follows 
from EguationsfTlandfTSlthat the accreted mass density required to 
produce the CIB fluctuations by redshift Zend is roughly 


Pace— 2x10®M© Mpc ® 


5^cib 


Aci 

0.2 


(s) 


0.072 nWm-2sr-i 

V10 A1-0 


(16) 


where e is the time-averaged efficiency, faed is a bolometric cor¬ 
rection factor i.e. the fraction of the total energy emitted that ends 
up being observed in the NIR (2-5pm). This requirement is illus¬ 
trated in Figure With maximal efficiency (e=0.4) we require 
a minimum of pacc>4x lO'^M© Mpc“® for producing the CIB 
fluctuations. If this is to be reached at high-z, any accretion ac¬ 
tivity by these black holes at later times would inevitably grow 
the net mass density locked in BHs further. Given the brief cos¬ 
mic time elapsed at high-z, this lower limit for the accreted mass 
is quite large in comparison with the mass density in SMBHs 
in the local universe z=0, estimated at (4—5) x 10®M©Mpc~® 


dShankar et alj2004lMarconi et alj2004lVika et alj2009lrl Other 
studies place additional limits on the mass de nsity at Z'^6 to 
<10^ M©Mpc~® based on the unresolved CXB dSalvaterra et al.l 
l2012h : and <10®M© Mpc~® based o n the integrated X-ray emis¬ 
sion of high- 2 ; sources jTreister et alj|201^ . Both these limits can 
however be avoided if the emerging SED is significantly altered, 
or if BH cores are heavily obscured or the accretion is radiatively 
inefficient. Relaxing the assumptions of constant e would like¬ 
wise alter these numbers. Nevertheless, for the CIB to be entirely 
from high -2 we require an abundant radiatively efficient popula¬ 
tion of black holes established very early on, with significantly 
slowed/inefficient growth over the ma jority of the remaining cos¬ 
mic time (see e.g. iTanaka et al.ll2012h . Whether the required pbh 
is in fact realistically attainable is a matter of debate and outside 
the scope of this paper. In what follows we do not need to make as¬ 
sumptions on the BH seed masses and growth mechanism but will 
instead assume that the required pace can be reached by 2=10 and 
explore the observational consequences. 



Figure 9. The requirements for accreted mass density and time averaged 
efficiency to produce the reference <5Tcib level. The lower limit corre¬ 
sponds to faed = 1 when all the energy emitted ends up being observed in the 
NIR. The requirements increase quickly as the spectmm becomes broader 
/sed<l- The solid lines correspond to 2eii(j=10 whereas the dotted lines 
show 2end=15,20 (lowei'/upper). The right veifical axis shows the accreted 
density as a fraction of the baryonic density. 



Energy keV 


4.1 Correlation with the Unresolved Cosmic X-ray 
Background 

The recently detected cross-correlation signal between the unre¬ 
solved CIB (3.6 and 4.5pm) and CXB (0.5-2keV) fluctuations 
was s uggested to originate f rom black holes in the early uni¬ 
verse dCapnelluti et alj I^13h . After acco unting for unr e solved 
AGN and X-ray binaries within galaxies, iHelgason et al.l d2014l) 
find that a tentative CIB x CXB signal (2.5a) remains at large- 
scales exhibiting the same shape as the CIB fluctuation power 
spectrum. From the definition of the spatial coherence C{q) = 
PcibxCxb/{PcibPcxb) and that {5 ~ q'^ P(q)/ 2tv one can 


^ Recent evidence su ggests this could be a s much as five times higher than 
previously estimated IComastri et al .120151) 


Figure 10. A simplistic broad-band SED in the rest-frame with two con¬ 
tributions 1) opticalAJV which is normalized to unity, and 2) a hard X-ray 
component which is a factor of 10“^ ,10”^, 10“^, 10“^ less energetic. For 
reference, the dotted curves show the emerging spectrum from the DCBH 
model of lYue et alJ b013bh also nomialized to unity in the UV/optical. 
These spectra are used as templates in Figure HTI 


write 

xCXB — VC SFcib SFcxb (17) 

where 5 TcibxCxb is the CIB x CXB cross rms-fluctuation. Taken 
at face value, the large scale average (3'-30') of this measure¬ 
ment gives 5 TcibxCxb~ 10 ^^ergs ^cm ^nW m ^sr ^.With 
the CIB fluctuations being roughly 0. 04nW m ^sr ^ 

(see Section|3 we must have v/C(5Tcxb— 10“®erg s“^cm“^sr“^ 
to satisfy the CIB x CXB correlation. From Equations [T3 and [T] the 
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Flux Fluctuations 




Figure 11. Left: The cumulative CIB (upper) and CXB (lower) from black holes at high -2 growing according to the BH growth in the u pper right panel. The 
pace at 2=10 has been chosen such that the measured 51^4.5^m=0.037di0.010 is reproduced (marked as K12, iKashlinskv et alJ 2012 h . The results here are 
mostly independen t of the growth history and the individual BH masses. The upper left panel shows the resulting -F 4 . 5 ;^m CIB where TeV upper limits from 
Meyeret aU j2^2h are also shown. The lower panels show the CXB assuming three SED tem plates for the X/Q=-2 , -3,-4 ( yellow, gree n, blue) and that from 
Yue et alj i2013bh (dashed curve) shown in Figure[To| Different upper limits for the CXB from lHickox & Markevitchl fcOOd) (HM06) and iMoretti et al.l fc012h 
(M13, converted to 0.5-2 keV flux) are shown as da shed lines with Icr unce rtainties as grey regions. The levels of the lai'ge scale CXB fluctuations are shown 
in the lower right panel with the measured levels of ICappelluti et al] j2013b (Cl3) shown for reference. Right: Accreted mass density in black holes growing 
at Eddington with e=0.1 reaching pbh=3x IO^MqMpc”^ by 2=10 (solid). The dashed line shows the lower limit for SFqiq to be pr oduced according 
Equation llbl in the extreme limit and e=0.4 . The m ass density locked in SM BHs at 2=0 is also shown for reference as grey regions jVika et al.l20Q9b 

(V09) and upper limits at 2~6 from ISalvaterra et*^ j2012h and iTreister et al.l i2013h in blue (upper, lower respectively). The lower panel show s the CIB x CXB 
fluctua tion s in units of erg s~ ^cm~^nW Notice how the JY/0=-2 case (yellow) is the only one to reach the measured levels of ICappelluti et al] 

<2Q13h and I Yue et al.l i2013bh (dashed) appears notably below this level. 


implied CXB flux is 


F’cxb=7x 10 ^^ergs ^cm ^deg 


(18) 


Since Fcxb is minimized for perfectly coherent sources, C=l, this 
value represents the amount of X-ray flux associated with the CIB 
sources whereas it is a lower limit for the total unresolved CXB. 
However, this flux can be accom modated by independent meas ure- 
ments of the unresolved CXB in (lHickox_&_Markeyitchl l2Q0^ but 
is only marginally consistent with Moretti et al.l ( 1201 2h who find 


only ^3x 10~^^ erg s~^ cm~^deg~^ to be still unresolvetE How- 
ever. lMoretti et al.l j2012l) do not include energies <1.5 keV in their 
fit for the CXB spectrum which could be the source of their lower 
derived 0.5-2 keV flux. On the other hand, the required Fx could 
itself be reduced if i) the large scale clustering of the underlying 
sources is stronger, say with Aci>0.2, or ii) if the CIB x CXB level 
at large scales is actually lower in reality, which is possible within 
the un certainties, bringing it into agreement with low -2 AGN and 
XRBs dHeleason et ^2014ll . 

To illustrate these requirements, we display four simplistic 
SEDs in Figure [TO] where we assume an optical/UV component 
and an X-ray component that is a factor of 10“^,10~^,10“®,10~^ 


® We have converted the l.SkeV flux of iMoretti et alj d2012h to 0.5-2keV 
using their spectral index r=0.1 


less energetic. We refer to these in terms of the X-ray to opti¬ 
cal ratio, X/0=—1,—2,—3 and —4 respectively where X/0= 
logj^gLx/io with Lx and Lo integrated in the ranges ( 1 -|- 
2 )[ 0 . 5 , 2 ]keV and [4.5/(l-|-2),0.091]^m respectively, with 2=10 
. These ratios are already significantly lower than locally observed 
AGN which exhibit X/O~0 on average but a suppressed X-ray 
component is possible for highly obscured BHs. The ionizing pho¬ 
tons are then absorbed by the thick gas and reprocessed into nebular 
emission in t he UV/optical (L y-a, free-free), enhancing the CIB 
contribution dYue et al.ll20I3bh . In this way, a larger fraction /sed 
of the bolometric luminosity ends up in the CIB, requiring less ac¬ 
cretion onto BHs (see Fig. H. We dis play the emerging SED from 
the highly obscured DCBH model of l Yue et H] d20^13bl) in Figure 
llOlfor comparison. The energy at longer IR/radio wavelengths, and 
shorter 7 -ray wavelengths are assumed to be negligible in compar¬ 
ison. 

Figure [TT] shows the cumulative CIB and CXB flux including 
fluctuations and CIB x CXB cross fluctuations resulting from these 
broad-band SEDs. Because the shape of the growth history does not 
matter much, we simply assume a scenario in which all black holes 
grow at the Eddington rate with e=0.1, reaching a fixed accreted 
mass density at 2 end =10 (=pio) 


Pace ( 2 ) = pro exp 


F l-e \ fio-fz 

\ e J fsdd 


(19) 


where fsdd—0.45 Gyr is the Eddington time scale. The results 
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are also not sensitive to the choice of Zend as long as we keep 
e and pio fixed. We have chosen pio=3xlO®M0 Mpc“® such 
that the SFcib is reproduced. We find that the indicative large- 
scale CIBxCXB signal can only be attained when X/0~—2. At 
the same time, the unresolved CXB implied by X/0>—2 is in 
slight tension (l-a) with the most stringent limits measured by 
Moretti et alj (I 2 OI 2 I 1 (see Fig. [TT] jottom left). The spectrum of 
YueetalJj2013bl) laving —4<X/0<—3 is shown as a dashed line 
in Figure [TT] and appears notably below the CIB x CXB require¬ 
ment. This means that the lYue et all (1201 3bll model actually falls 
short of the X-ray energy implied by the coherent CXB levels but 
is nevertheless consistent with the CIB levels. 

It is important to emphasize that the sources producing the 
CIB and CXB are not necessarily the same physical emitters and 
are in general a mixture of stars and BHs sharing the same large 
scale structures. However, if they are predominantly BHs emit¬ 
ting in both optical and X-rays they will be accompanied by ioniz¬ 
ing UV emission. For AGN spectra, this contribution is estimated 
by connecting the optical and X-ray part of the spectrum with a 
power law slope, shown as dashed lines in Figure [Toj If we adopt 
the measured amplitudes of the CIB and CXB fluctuations we ob¬ 
tain a slope ranging from dlog(Fb /Fx)d\ogE=—0.5 to —0.7 and 
the resulting ionization rate exceeds the recombination rate in all 
cases such that freionrecomb- TMs meaus that in the presence 
of such production rate of ionizing photons at z>10, the Universe 
will reionize in the matter of ~tens of Myr. For the CIB x CXB 
signal to arise from the same type of sources therefore requires the 
suppression of UV photons esca ping into the IGM ( . feBc<l) e.g. the 
black holes to be ultra-obscured dRicotti et al.l200a : IComastri et al.l 
I2OI5I) . This would be easier with unconventional BHs such as 
DCBHs rather than X-ray binaries from stellar remnants which 
are much less obscured and are limit ed in growth by the mass and 
lifetime of the companion star (e.g. I Yue et al.ll2^ : iTanaka feuil 
I 2 OI 4 I : IPacucci & Ferrarj|2015h . Also, for an abundant population 
of small BHs, a proportionally greater fraction of the net mass goes 
into the seed of the gravitational collapse and does not contribute 
to Pace. 


5 SUMMARY AND DISCUSSION 

In this paper, we assessed the possibility that the established clus¬ 
tering signal in the source-subtracted CIB fluctuations at near-IR 
originates in the early universe. We consider emission from star 
formation and black hole accretion during the first 500 Myrs and 
establish the physical requirements and consequences to reach the 
measured CIB values. Our results can be broadly summarized as 
follows: 

• The contribution from high-z galaxies, with star formation 
efficiency set to /* =0.005 to match the LF from the deep HST 
Legacy surveys, cannot reproduce the levels of unresolved fluc¬ 
tuations in the near-IR for any reasonable extrapolation, pro¬ 
ducing merely ~0.01—0.05nWm~^sr~^ flux at z>8. In fact, 
faint galaxies at intermediate-z make a larger contribution but are 
themselves insufficien t to account for the entire clustering signal 
dHelgason et al.ll2012h. This conclusion is in agreement with previ¬ 
ous stu dies dFemandez etal]|2012l : ICoorav et al.ll20123 : lYue et alj 
l2013ah . 

• A simple extrapolation of the faint-end of the LF suggests that 
high-z galaxies will become significant once sources have been 
removed down to mAB^32. Using the CIB to probe these popula- 
tions will therefore require JWST dKashlinskv et alj|2015tl and/or 


cross -correlation techniques dAtrio-Barandela & Kashlinskvl 

|2014) . 

• For more energetically efficient modes of star formation, we 
find that heavier IMFs still require high star formation efficien¬ 
cies /*>0.1 which must either take place in small halos or end 
before z~12. The ionizing photons associated with such a popu¬ 
lation would have to be suppressed to avoid rapid reionization and 
to enhance the nebular emission contributing to the CIB. In other 
words, it would require very low escape fractions , /esc<0.1%, in 
low-m ass systems towards high redshifts (see also [Fernandez et al.l 
bold). This is contrary to what is found in theoretic al studies 
IFerrara & Loebl[2013l : iKuhlen & Faucher-Giguerj|2012l) . 

• The cumulative CIB contribution of the fossils of the first 
galaxies that contain old stellar populations radiating throughout 
cosmic history is comparable to that of their high-z progenitors. 
Since many such systems will merge or become satellites of larger 
galaxies, their net contribution to the source-subtracted CIB is 
smaller and can be safely neglected as long as the contribution of 
the high-z counterparts is also low. 

• For the CIB fluctuations to be produced by accreting black 
holes, one requires vigorous accretion rates in the early universe 
reaching pacc>lO®M 0 Mpc“® by z=10. This quantitative limit 
assumes a constant efficiency but is mostly independent of the ac¬ 
cretion rate history and BH seed masses. The accretion must be 
very radiatively efficient at early times and then drop to inefficient 
growth throughout most of the remaining cosmic time with the re¬ 
sulting populat ion largely unacc ounted for in the SMBH census in 
local universe dVika et al.ll2009h . In addition, this population must 
be extremely gas-obscured in order to i) avoid reionizing the uni¬ 
verse too rapidly, ii) to suppress the associated soft X-ray back¬ 
ground, iii) to enhance the fract ion of energy emi tted in the UV 
as reprocessed nebular emission. lYue et all ( 12013bh have designed 
a novel model of DCBH formation that is broadly consistent with 
these requirements. 

• For the expected level of clustering of sources (both at low 
and high-z), the minimum soft CXB flux implied by the measured 
CIBxCXB cross-correlation, Fx—7x10“^® erg s“^cm“^deg“^, 
is only marginally c onsistent with the most stringent CXB limits 
dMoretti et al .11201^ . At face value this means that either i) that 
high-z BHs require the entire available unresolved CXB, ii) the 
clustering of the underlying sources is stronger (e.g. more biased) 
with Aci>0.2 at sub-degree scales, or iii) the high current un¬ 
certainty of the CIBxCXB level at large scales c ould eventually 
bring it into agreement with low-z AGN and XRBs dHelgason et al.l 
|2014 . 

In this paper, we explored the physical requirements for stellar 
emission and accretion separately. It is however entirely possible, 
and in fact plausible, that the CIB and CXB are not produced by 
the same physical sources. The most intuitive scenario would have 
the CIB produced by stars, whereas accreting black holes sharing 
the same structures are responsible for the CIB x CXB correlation. 
This would be a more natural way to obtain the Fir/F x~ 100 ratio, 
which is unconventional for BHs only. The requirements for BHs 
to produce the CIB x CXB component only, either stellar remnants 
or DCBHs, are much less energetically demanding than for the 
single-population explanation of the CIB and CIBxCXB fluctua¬ 
tions. However, whereas rapid reionization is more easily avoided, 
one still requires the entire unresolved CXB available and is then 
left with the problem of explaining the CIB, the subject of this pa¬ 
per, by other means. 
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Our study of stellar sources isolates a narrow comer of pa¬ 
rameters required for explaining the properties of the measured 
CIB fluctuations in terms of physics and evolution at high-2. The 
predominant contributions would then have to come from very 
massive stars (>100Mq), and accreting black holes, with non- 
instantaneous reionization in a highly clumped IGM. The area cov¬ 
ered by this narrow corner may be widened by a factor of a few by 
relaxing assumptions of constancy of parameters such as /* or e. 
Apart from nucleosynthesis and accretion, an additional emission 
component at high 2 may come from stellar collisions in the dense 
stellar system phase which ma y result from stellar dynam ical evolu¬ 
tion of the early systems (e.g. iKashlinskv & Re^l 19831) . The spe¬ 
cific weight of this component is hard to quantify in a robust model- 
independent manner, but it may result in additio nal energy releases 
and u ltimately formation of supermassive BHs JBegelman & Reed 
Il978h . It is important to emphasize that narrow allowed parameter- 
space is predicated on 1) the assumption of accurate estimate of the 
power contained in the primordial density spectrum at very small 
scales (comoving kpc), and 2) the assumption of its time-invariance 
after the first sources form. The first drives the collapse of the first 
halos with tiny modifications in the amount of power leading to 
exponentially different fcoi{z) and, consequently, different, poten¬ 
tially larger CIB. The second possibility may arise from the in¬ 
evitable and highly complex interaction between the condition of 
the density field which form s early collapsing objects and the en - 
ergy releasing sources (e.g. lOstriker & Cowielll98ll : lRee ^ Il985h . 
In this context, recent studies have found that high-2 galaxies may 
have had to form surprisingly early to account for the weakening 
of the expo nential decline at t he bri ght end of the LF with increas¬ 
ing redshift lpinkelstein et al.l ( l2015h . This could be achieved by an 
unusually early collapse with respect to standard clustering theory, 
or an increased star formation efficiency making t he stellar-to-halo 
mass ratio appear to rise with increasing redshift dSteinhardt et al.l 
l2015h . In order to have a significant impact on the CIB however, 
these effect must also be present in lower mass systems which con¬ 
tribute the bulk of the CIB. 

All this argues for empirically measuring the CIB fluctuations 
with required fidelity, establishing its nature and then constructing 
and identifying the implied physics models for the sources pro¬ 
ducing the signal with all its properties. This has been demon¬ 
strated in multi-wavelength cross-correlation studies that provide 
further constrain ts for theoretical models of both low-2 and high- 
2 sources (e.g. ICappelluti et alj|2013l : iThacker et ai]|2014h . Ad¬ 
ditional information on the physics and evolution at early times 
would co me by looking for correlations betw een the CIB and 21 
cm maps iFemandez et al.ll20l4 lMaoll2014h . Euclid’s large sky 
coverage and wavelengths make it particularly suitable for us¬ 
ing its to-be-measured CIB fluctuations to isolate emissions from 
the epoch of re-ionization (Kashlinsky et al. 2015, submitted) 
and identifying the condition of the IGM by cross-correlating Eu- 
clid’s measured source-subtracted a ll-sky CIB and CMB maps 
dAtrio-Barandela & Kashlinsky! l2014h . A specific experiment has 
been recently proposed with the /IFST/NIRCAM to probe the 
Lyman-break of the CIB compone nt and to tomographical ly recon¬ 
struct the emissions out to 2>30 dKashlinskv et alJl20I5tl : the Ly¬ 
man tomography method proposed there was shown to already lead 
to interesting limits at 2>30 from the available Spitzer CIB maps 
at 3.6 and 4.5p.m. The nature of the CIB fluctuations should be re¬ 
solvable in the coming decade and the emissions produced during 
the epoch of reionization identified from its properties, shedding 
light on the physics and evolution at that time. 
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